Can buckling account for the features seen in graphite's Raman spectra? 
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Raman scattering data were collected on graphite monochromator. Spectra were interpreted in 
terms of the space group PQzmc, a subgroup of space group P63 frame. The latter has commonly 
been used for the interpretation of Raman scattering data. Space group PQ^mc corresponds to 
the buckling of graphene sheets and is consistent with many spectral features. Both the first and 
second order scattering were considered. Many first order results (most notably the assignments 
of the band at 1350 and the peak at 1620 cm -1 ) were found to agree with previous observations 
[Y. Kawashima and G. Katagiri, Phys. Rev. B 66, 104109 (2002)] carried out on highly oriented 
pyrolytic graphite samples. To check the consistency of the model, symmetry analysis was applied 
to the second order spectra. Also a simple test for buckling model was done. 
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I. INTRODUCTION 

Despite its nominally simple crystal structure the 
interpretation of the Raman and infrared spectra of 
graphite has been an active and controversial topic. The 
crystal symmetry of graphite is reported to be either 
P6^/mmc (No. 194) or PQ^mc (No. 186)^. A more rare 
form, with rhombohedral space group symmetry R3m 
(No. 166) was also reported in ref. OJ The symmetries of 
these different forms increase in order Rim — > P&^mc — > 
P63 /mmc: these symmetry groups share group-subgroup 
relationships. It is well known that, despite its virtually 
simple structure, a large number of features which seem 
to contradict the space group symmetry PQ^/mmc, have 
been observed in Raman spectral. Partially the prob- 
lem has been that rather different types of 'graphites' 
were studied, as far as crystallite size and defects are 
considered. Despite that Raman spectra were commonly 
interpreted in terms of space group P63 / 'mmc. 

Still, the puzzling features of the Raman spectra of 
graphite are (i) the peak(s) at around 1355 cm -1 (of- 
ten labelled by D), (ii) peak at around 1620 cm -1 (of- 
ten labelled by D') (iii) the first overtones of the afore- 
mentioned modes (often labelled by 2D and 2D', respec- 
tively). If PQ^/mmc symmetry is assigned to graphite, 
the Brillouin zone centre modes transform as the irre- 
ducible representation (IRREP) 2A 2u © 2B 2g © 2E iu © 
2E 2g , where A 2u and E\ u representations are infrared 
active, B 2g is silent and E 2g is the only Raman active 
representation. Thus, only two peaks, observable at xx, 
yy, xy and yx geometries, are expected (i.e., E 2g modes 
can only be observed if the polarization vectors of both 
the incoming and scattered light are perpendicular to the 
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hexagonal c axis). They were reported to appear at 42 
and 1580 cm -1 (the latter is commonly labelled as 'G' 
mode) 3 . However, the number of observed peak is gen- 
erally greater. Now, it is interesting to note that in the 
case of PQ^mc symmetry the Brillouin zone centre modes 
transform as the IRREP 2Ai®2Bi®2Ei®2E 2 , where A x 
and Ei symmetry modes are Raman and infrared active, 
Bi mode is silent and E 2 mode is Raman active. 

The first attempt to explain the D band were based on 
disorder, which was interpreted to cause an appearance 
of a totally symmetric modai As is discussed below, 
this mode(s) is not totally symmetric and its origin is 
still not clear. Recent explanations are based on dou- 
ble resonance ideaS&i. The merits of the double reso- 
nance model are its ability to take the laser beam wave- 
length dependent positions of D and 2D bands into ac- 
count. The first attempt to explain the D band were 
based on an idea that, due to the disorder, a Raman in- 
active mode became active^. The intensity of this mode 
was reported to be dependent on particle sizei In ad- 
dition, the frequency of this mode was reported to in- 
crease with increasing laser light frequency^. The corre- 
sponding changes occurred in the first overtone of the D 
band. To explain these features a double resonance Ra- 
man scattering model was develope d 2 ! 6 ' 7 . Several weak 
peaks were observed in Raman spectra, the most no- 
table at 867 cm -1 , which was reported to correspond to 
the out-of-plane vibration*!^. This model assumes that 
the D and D' bands are transverse optic Brillouin zone 
boundary [K point) and zone centre near Brillouin zone 
centre modes. The double resonance model further ne- 
cessitates that the atom displacements occur in graphene 
sheets (in-plane vibrations) 6 . Occasionally, double reso- 
nance theory has been merited to take the slight differ- 
ence in the anti-Stokes and Stokes positions into account. 
In this context once should bear in mind that, in order 
to observe the high frequency overtones in anti-Stokes 
spectra, exceptionally high laser beam power densities 
were used in these studies (the highest powder densities 
were approximately 5000 times larger than the one used 
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in this study Although this type of considerations are 
interesting, we note that, in the case of the anti-Stokes 
spectra, the reported differences between the Stokes and 
corresponding anti-Stokes frequencies were, for different 
modes, positive, almost zero and negative, as was pointed 
out in ref. la 

The purpose of this work was to reinvestigate a high 
quality graphite sample by Raman scattering technique. 
The first goal was to check if the buckling model can ex- 
plain the first and second order spectra. The second goal 
was to consider to which extent the proposed double res- 
onance model is consistent with the observed symmetry 
properties of Dand D' bands. 



II. EXPERIMENTAL 

A pyrolytic graphite monochromator (Advanced Ce- 
ramics) sample was used in this study. Raman mea- 
surements were performed using a Jobin-Yvon T64000 
spectrometer consisting of a double monochromator cou- 
pled to the third monochromator stage with 1800 grooves 
per millimeter grating (double substractive mode). Prior 
to measurements, spectrometer was calibrated with Ne 
lamp. Acquisition time was adjusted to have a sufficient 
signal-to-noise-ratio. CCD detector was used to count 
photons. Backscattering measurements were carried out 
under the microscope. Also 90-degree scattering angle 
experiments (the angle between the incoming and scat- 
tered light was 90 degrees) were carried out. Raman spec- 
tra were excited using an argon ion laser. In the case of 
backscattering experiments, the laser beam power on the 
sample surface was 200 /iW and the diameter of the laser 
beam spot was approximately 2 /zm (the spot diameter 
was approximately 100 /im in 90-degree scattering angle 
measurements, and the power was 50 mW, so that the 
powder density was smaller than in the case of backscat- 
tering measurements). In a backscattering geometry the 
wavelength was 514.532 nm and in a 90-degree scattering 
angle measurements it was 487.986 nm. 

III. RESULTS AND DISCUSSION 

Selected regions of Raman spectra collected on aft- 
basal plane with an polarizer set parallel (labelled as 
Z(XX)Z) and perpendicular (labelled as Z(XY)Z) with 
respect to the incoming polarization direction are shown 
in Fig. G](a), (b) and (c). Particular attention was paid 
on the D, D', G, 2D, and 2D' modes. Here Z and Z 
refer to the incoming and scattered light propagation di- 
rection. Although the light propagation direction can be 
rather accurately determined at backscattering geome- 
tries, incoming light polarization direction with respect 
to hexagonal a and b axes was not known. The corre- 
sponding regions, obtained through right-angle geome- 
try, are shown in Fig. H(d), (e) and (f). The intensity of 
G mode at around 1580 cm -1 was the same for XX and 



XY geometries (see Fig. 01 ( a ))> consistently with the 
idea that it belongs to the E 2g symmetry (space group 
PQ^/mmc) or E% symmetry (space group PQ^mc). In 
contrast to the earlier reports (see, e.g., ref. 3), we could 
not see any other mode below 1580 cm -1 . Fig. ^ (b) 
shows the broad band centered at around 2750 cm -1 . 
We note that this band is consisted of several peaks and 
the intensity differs from zero in XX, XY, ZZ and ZY 
geometries. Now, it is of interest to consider the symme- 
tries of the first overtones (second order scattering) in the 
case of the fundamentals belonging to the point group 
symmetries 6/mmm and 6mm. Only two-dimensional 
IRREPs need to be considered, since the first overtone 
of a one dimensional IRREP is always totally symmet- 
ric and can be distinguished by depolarization measure- 
ments. Now, the first overtone of the E\ u , E± g , E 2u and 
E 2g modes transform as A\ g © E 2g and the first over- 
tone of the E\ mode transforms as A\ ®E 2 . Thus, if this 
band is due to the second order scattering of the Brillouin 
zone centre mode, it must belong to E\ u (space group 
P63 / mmc) or to E\ symmetry (space group P6^mc) . To 
check if either assignment is correct two other experimen- 
tal geometries were used. First one is the geometry where 
light propagates perpendicularly to the hexagonal c-axis, 
Fig. n (d)-(f). Consistently, G mode was significantly 
weaker (if the polarization direction of the incoming or 
scattered light is strictly parallel to the hexagonal axis, 
the G mode should have zero intensity. Deviations might 
be because of the folding of graphene sheets). The ap- 
pearance of the mode at around 1355 cm -1 is consistent 
with the space group P63TOC. Basically, if this mode be- 
longs to the Ei symmetry, it should only be observed at 
the XZ and YZ geometries. Thus, this mode should not 
be observable at ZZ geometry, in constrast to the exper- 
imental observations. Although partially this might be 
due to the aforementioned folding we decided to carry 
out 90-degree measurements. 

If the buckling model alone is sufficient for explaining 
the appearance of D and D' bands (instead of the dis- 
order, which activates non-Brillouin zone centre modes), 
one should see them at 90-degree measurements carried 
out on large graphene sheets (so that the effect of dis- 
order as a dominant factor can be eliminated). Now, 
the phonon propagates parallel to the graphene sheet, by 
momentum conservation rule. The G mode was observed 
(Fig. I2J), since both the incoming and right-angle scat- 
tered light have a polarization component parallel to the 
ab plane. The more interesting point is that now also the 
modes belonging to the E\ symmetry modes should be 
observed. Thus, this experiment served as a test to clar- 
ify whether D and/or D' modes are due to the Brillouin 
zone centre modes or if they correspond to the disorder 
activated modes. As Fig. [21 reveals, D and D' bands 
were not observed, but 2Dwas observed. Rather com- 
mon opinion is that D mode is related to the Brillouin 
zone K point mode (sometimes assigned to the totally 
symmetric mode) , whereas D' mode has sometimes been 
considered to belong to B\ symmetry with wavevector 
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FIG. 1: Raman spectra collected on ab plane at Z(XX)Z and Z(XY)Z geometries: (a) spectral region corresponding to the 
D, G and D' modes, (b) Spectral region corresponding to the 2D mode(s) and (c) 2D' modes. The corresponding Raman 
spectra collected on edge plane at X(ZZ)X and X(YZ)X geometries are shown in panels (d), (e) and (f). Spectra shown in 
red correspond to the setting where polarizer and polarizator were parallel to each other, whereas blue lines give the spectra 
collected on with polarizer perpendicular to the incoming polarization direction. 
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FIG. 2: Raman spectra collected on 90-degree scattering 
angle geometry, where phonon propagates along the graphene 
sheet, (a) spectral region between which should reveal also 
the D and D' modes, (b) Confirmation that 2D are observed 
at this geometry, too. 



close to zero (Brillouin zone centre), see for example ref. 
2. Now, the problem is that these assignments do not 



comfort with the polarization measurements. The first 
overtone of these modes should also be totally symmetric, 
in contradiction with the experiments. Now, the selec- 
tion rules for phonons which can participate on double 
resonance scattering state that only totally symmetric 
modes can resonantly couple electrons within the same 
non-degenerate bands 2 . The second possibility is that 
the phonon involved in the process should belong to the 
Bi symmetry. 



IV. CONCLUSIONS 

The Raman spectra collected on graphite monochro- 
mator were interpreted in terms of two space group sym- 
metries. The only observed first order lines were G mode, 
and D and D' modes. It was concluded that although 
the weak, additional lines reported in literature probably 
do correspond to the buckling of graphene sheets, the 
D and D' modes are probably not due to the buckling. 
The latter feature was confirmed by 90-degree angle ex- 
periments. However, the situation is different for smaller 
crystal size materials. It was also confirmed that D and 
D' modes do not belong to the totally symmetric repre- 
sentation. 
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